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ABSTRACT 

The Milky Way harbors giant H II regions which may be powered by star com- 
plexes more luminous than any Galactic OB association known. Being across the 
disk of the Galaxy, however, these brightest associations are severely extinguished 
and confused. We present a search for one such association toward the most lu- 
minous H II region in the recent catalog by Murray and Rahman, which, at ~9.7 
kpc, has recombination rate of ~ 7 x 10 51 sec -1 . Prior searches have identified 
only small scale clustering around the rim of this shell-like region, but the pri- 
mary association has not previously been identified. We apply a near-infrared 
color selection and find an overdensity of point sources toward its southern central 
part. The colors and magnitudes of these excess sources are consistent with O- 
and early B-type stars at extinctions 0.96 < A K < 1.2, and they are sufficiently 
numerous (406 ± 102 after subtraction of field sources) to ionize the surrounding 
H II region, making this a candidate for the most luminous OB association in the 
Galaxy. We reject an alternate theory, in which the apparent excess is caused by 
localized extinction, as inconsistent with source demographics. 

Subject headings: open clusters and associations: individual — stars: massive - 
stars: formation — infrared: stars 



1. Introduction 



Very massive OB associations are objects of intense interest. Given that the association 
birthrate (J\ r \) is observe d to vary with association mass (M c \) as dN c \/d\riM c \ cx 1/M C \ 
jMcKee fc Williamsl fl9~9~7T ) . or possibly slightly flatter, a large portion of each galaxy's star 



formation occurs within its most massive OB associations. Giant OB associations are most 



-2- 



capable of disrupting their gaseous environment, from their natal molecular clouds to the en- 
tire galactic neighborhood, and they inflate superbubbles which erupt from the disk, feeding 
a galactic halo or fountain. They sculpt and illuminate the giant H II regions seen in distant 
galaxies and are a notable feature of starburst evolution. Being the most extreme examples 
of intense star formation in the current universe, they are laboratories for the physics of star 
cluster formation such as environmental influences on the initial mass function, formation of 
very massive stars, and dynamical evolution with rapid stellar evolution. 

Although several of the brightest know n OB associations within the Milky Way approach 
'super star cluster' status (M > 1O 4 M , iPortegies Zwart et al.ll201Q[ ). it is possible that 
even more luminous associations have escaped detection because of severe extinction and 
confusion in the Galactic plane. For instance, the Galactic OB associations with the largest 
known ionizing output, S, are those powering the NGC 3603 H II, Arches, Quintuplet 
and Galactic Center reg ions (logS'/s -1 ~ 51.5,51.0,50.9, and 50.5, respectively: iFigerl 12008 



Conti fc Crowtherll2004j ); however, the upper limit o f the Milky Way associati on distribution 
is estimated to be even brighter (logS^/s -1 ~ 51.7; lMcKee fc Williams! Il997| ). 



Recently, iRahman fc Murray! (120101 ) identified 40 star-forming complexes within the 
13 most lumi nous (S > 3 x 10 51 s -1 ) fr ee-free emission regions observed by WMAP and 
catalogued by iMurray &: Rahman! (120101 ) . Many of these complexes are potential hosts to 
embedded super star clusters. Of particular interest is the most luminous unconfused star- 
forming complex located at (/, b) = (298.4°, —0.4°) with a kinematic distance of 9.7 kpc. 
With an output of logS/s -1 ~ 51.8, the stars powering this region may represent the Milky 
Way's most luminous OB association. (We refer to this region and its cluster as G298 or, 
inspired by its 8/im appearance, as the 'Dragonfish Nebula': Figured!) The region appears 
in 8 /im as a closed bubble surrounded by a prom inent shell. Previo us searches using near- 
i nfrared (NIR) poin t source catalogs from 2 MASS JDutra et alil2003h and Spitzer GLIMPSE 
(IMercer et al.ll2005l ) have concentrated on the brightest regions of free- free emission around 
the shell and relied on integrated source counts. Although some clustering was discovered 
in these studies (Figure [2]), the central ionizing association is yet undiscovered. 

In this Letter, we identify a candidate for the central OB association powering the 
Dragonfish, primarily by applying a color selection to the 2MASS catalogue. Our candidate 
is potentially the most luminous (and thus massive) OB association in the present Milky 
Way. 
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Fig. 1. — G298 or the Dr agonfish region seen in the Spitzer GLIMPSE 3.6 (blue), 4.5 (green) 
and 8.0 (red) /im bands ( jBenjamin et al.l 120031 ) . 
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2. Feasibility and Methods 



Only very luminous stars of deeply embedded OB associations are visible due to distance 
and extinc tion. Whether the a ssociation can be identified using the 2MASS point source 
catalogue (ISkrutskie et al.ll2006l ) depends on the OB association radius and field star density 



Radii in the range 10-30 pc ar e likely since m ost associations, irrespective of mass, are 
disrupted immediately after birth ([Roberts! I1957D . OB association birthsites in the Milky 
Way have a limited range of mass column densities, S ~ 0.3 gem -2 ± 0.5 dex ( jFall et al. 
2010l ). Turbulent linewidths in these birthsites are a = 5a 1 / 2 [(0.5/£:)(X/.3gcm _2 )M5] 1 / 4 
kms -1 , for associations of mass 10 5 M 5 M which form with efficiency e from a region with 
turbulent virial parameter a. Crossing and formation timescales are both shorter than 
the ionizing lifetime and the terminal expansion speed of an OB association is likely to 
be ~ a. For an association of ~ 10 5 M Q , an age between 2-4 Myr, and an expansion speed 
between 4-6 kms -1 , we therefore expect an expanding population currently 8-24 pc in radius, 
possibly surrounding a dense core which survived the disruption event. For a physical radius 
R = 10(i2/10pc) pc with the ionizing luminosity of the Dragonfish as measured from WMAP, 
the association source density is 10(i?/10pc) -2 arcmin -2 . We measure the 2MASS source 
population density in the Galactic plane to be 14 — 17 arcmin -2 . For R < 8 pc the association 
dominates raw 2MASS source counts, whereas for R > 10 pc it is a minor variation of the 
2MASS population which is difficult to detect using raw source counts. 

The NIR color information provides an additional avenue to search for the embedded 
associations. We make use of this in the search for the embedded association within the 
Dragonfish. Ionization from a coeval stellar population drops sharply after 4 Myr as its most 
massive and powerful stars turn off the main sequence, so it is unlikely, though not impossible, 
that G298's i onizing population is any older. The ionizing population is sufficiently rich 



(S > 10 50 s l : iKennicutt et al.lll989l ) to fully sample the stellar initial mass function (IMF); 



its very massive stars should still be present and visible. 

The association is likely highly extinguished, even if it has cleared its natal cloud. At 
a distance of 9.7 k pc the fiducial distance-ext inction relationship within the Galactic disk, 
Ay I D ~ 1.6kpc -1 ( Binney fc Merrifield 1998), pre dicts Ay ~ 16; t his corresponds to a K s - 



band extinction Ax s — 1-0 ( INishiyama et a 



l2008h. We adopt t he iNishivama etaD J2009h 



extinction law, which is steeper than that of ICardelli et al.l ( 119891 ). We rely primarily on the 
2MASS point source catalog for our analysis, but use three additional data sets to check for 
consistency: the USNO-B catalog of visible magnitudes, the Spitzer GLIMPSE catalog at 
3.6 and 4.5 /im, and deeper J and H band photometry of a small section of the candidate. 



For our search, we target main- sequence, massive stars, at a distance of ~ 9.7 kpc and 
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an extinction of around Ak — 1- Further, association members visible in all three bands 
have spectral types of roughly Bl or earlier, so their intrinsic emission in the NIR bands 
(A ~ 1.0 -2.4 mm) is entirely in the Rayleigh- Jeans limit ((J -H,H - K) = (-0.11, -0.10); 



Martins fc Plezll2006[ ) . With the predicted extinction, NIR colors should be ( J— H, H— K) 



(1.1,0.5). Association sources should lie in a tight grou p along the red dening vector from 



this common initial color. If the upper IMF follows the ISalpeterl (119551 ) slope to an upper 



limit of 120 M , we expect ~ 400 x 10~ a34( ^~ L0) (D/9.7 kpc) L16 members with confident 
NIR colors ([J,H,K S ] < [15.9, 15.0, 14.3], for which the signal-to-noise ratio exceeds five in 
all bands). Note that if the association is older than 4Myr or is not coeval, it must be even 
more numerous because more stars are required to make up for the missing ionization from 
the top of the main sequence. 

The first steps of our analysis are to select sources with the color of an early-type 
star behind a range of extinctions bracketing the expected value, to test for a statistically 
significant overdensity within the confines of the 8 /im bubble, and then to vary the adopted 
extinction range in order to optimize the statistical significance of the overdensity. In practice 
we determine source density using the angular offset 6 N between a source and its Nth. 
neighbour, as = N/ (n0%), and use N ~ 20 — 100 to reduce Poisson noise at some expense 
in resolution. Once the angular scale of the putative association is determined, we assess its 
significance by comparing source density within this region to the mean and standard 
deviation us of those within independent, identically-sized regions in a 2-degree field centered 
on, but excluding, the 8//m bubble. Statistical significance is measured by (£* — fX)/ox;- 

Once a candidate is identified, we check that the same procedure applied to regions 
outside the 8 /im bubble yields no candidates of similar significance. We then examine the 
color-color and color-magnitude diagrams (CCDs and CMDs) of similarly-sized regions to 
check that an apparent overdensity is not caused by features of the distribution of stars and 
extinction. 



3. Results 

3.1. Candidate OB Association 

Based on the expected source colors from §(2j we apply a NIR color cut of 1.0 < J — H < 
1.4 and 0.44 < H — K s < 0.62, corresponding to colors of hot stars with 0.9 < Ak < 1-2. 
Within the color cut, the 2-degree diameter field encompassing the Dragonfish has a mean 
stellar density = 1.4 arcmin -2 , which we refer to as the field star density. The standard 
deviation of X* averaged over association-sized regions is cr^ = 0.30 arcmin -2 . We find 
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a significant overdensity of stars in the Dragonfish Nebula at (I, b) = (298.55°, —0.72°), 
located inside the shell of the star forming complex (Figure |2]). The 1—a stellar overdensity 
contour (E* = + a^) has semimajor and semiminor axes (11', 10'), corresponding to 
(31, 28)pc; however, the association may extend beyond the observational limit. The peak 
of the overdensity is £* = 3.7 arcmin -2 , 8-cr above the background level on cluster-sized 
regions, adopting the 100th nearest neighbor. The region contains 897 sources within the 
adopted boundaries. Subtracting the contribution due to background sources, the candidate 
association is composed of 406 ± 102 sources, consistent with the value determined based on 
the chosen IMF and measured luminosity (§[2]). Using the 1—a boundaries, the mean density 
of the association is S* + 4.0<T£. We refer to the candidate as the "Dragonfish Association". 

The association is strongly concentrated within the central 4.1', but contains an extended 
asymmetric envelope. The central, dense, symmetrical region has a stellar density which 
depends as 9~ a with a = 0.44 ± 0.03 for 2' < 9 < 3', where the value of a depends on the 
assumed density of field sources. At larger radii the azimuthally averaged density drops more 
rapidly, but is dominated by the asymmetries visible in Figure 1. Half the excess sources are 
found within 4.7', or 13 (D/9.7kpc) pc, s imilar to the half-light radii of other Galactic OB 



associations (IPortegies Zwart et al. 



20101). 



Assuming half the projected mass falls within this radius, and assigning a total mass 
of 10 5 Mq we infer that the association's central regions are older than their virial crossing 
times, whereas its outskirts are not (T rr = 22 M yr). It therefore seems likely to be unbound 



on the whole ( IGieles fc Portegies Zwartl 1201 if ) but to contain a bound core; however this 



bears further investigation. 

Using the USNO and Spitzer GLIMPSE catalogues, we are able to confirm that the 
stars identified in the Dragonfish Association have magnitudes and colors consistent with 
OB stars with the stated extinction, but not further able to distinguish the association 
population from the background. In the GLIMPSE wavelengths (> 3 /mi), this is because 
nearly all stars have identical, Ray leigh- Jeans colors, whereas with the USNO wavelengths 
(< 1 yum), the spectral energy distribution is strongly dominated by the extinction. Thus 
template colors of most stellar types can be fit to the additional photometry by varying the 
assumed distance or extinction of the star. 



3.2. Color-Color and Color-Magnitude Diagrams 



The NIR color information provides an additional avenue to determine important pa- 
rameters of the candidate association, including the statistical significance of the overdensity, 
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Fig. 2. — Top: The 8 /im GLIMPSE image of the Dragonfish Nebula. The dotted ellipse 
indicates the location of the WMAP source, and the solid ellipse indicates the location of 
the star-forming complex. The cyan triangles indicate the location of known H II regions. 
Clockwise from top: G298.56-0.ll, G298.2-0.3, G298. 19-0.78, G298.9-0.4, and G298.837- 
0.347. The red squ ares and diamonds indica te the regions of stellar clustering identified by 



Dutra et all (12003( 1 and Mercer et al.l (12005( 1 . Bottom Left: The individual 2MASS sources 



within the field surrounding the Dragonfish Nebula passing the color cut. Bottom Right: 
The 2MASS stellar densities of the field surrounding the Dragonfish Nebula at the color cut 
corresponding to O-stars with extinctions between 0.9 < Ak < 1.2. The color bar indicates 
the stellar density in stars per square arcminute. 
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the possible background contributors, and the color variation of all point sources within the 
association boundaries. Figure |3] presents the CCD and CMD of the point sources within 
the 2 degree field centered on the Dragonfish Nebula, showing the distribution of all 2M ASS 
point sources, as well as the relative location of the color cut. We indicate in red the point 
sources that fall within the association boundaries. We indicate in blue all stars within the 
field that meet the association's color cut. The position of the main sequence is indicated 
at a distance of 9.7 kpc and e xtinguished by Ak = 1 .0. The reference magnitudes for the 



O sp ectral type are taken from lMartins fc Plea (120061 ). and for the remaining spectral types 



from iPicklesi (119981 ) . The path taken by an 05V star if placed at varying distances, assuming 
the average extinction to distance ratio, is indicated by the red line. The path taken by any 
given reference star is just a translation of this path on the diagram. The color cut is suffi- 
ciently large to permit spectral classes of A and earlier at the given reddening. However, the 
sensitivity of 2MASS limits stars at this distance and extinction to be no later than B1V. 

In Figure H] we bin the CCD and CMD to investigate the overabundance in the region 
within the Dragonfish Association as compared to the surrounding regions. These figures 
indicate the location and prominance of the cluster overdensity as a function of both color and 
magnitude. For the CCD, we bin the sources in a 35 2 bin grid spanning —0.3 < J — H < 2.5 
and —0.2 < H — K < 1.1. For the CMD, we similarly use a 35 2 bin grid with 5 < J — H < 
17.0, and identical J — H limits. We grid the point sources that fall within the association 
boundaries and subtract the average bin value from the surrounding field. In both the CCD 
and CMD, an overabundance of point sources appears in the range 1.0 < J — H < 1.5. We 
note an underdensity in the CCD at (J — H, H — K) = (0.8,0.3). The total number of 
sources "missing" in the underdensity is approximately 200. The underdensity is associated 
with a galactic structure feature in the surrounding field which we discuss in depth in a 
forthcoming paper. 

Given that the colors of these underdense sources are separated from those of the pu- 
tative OB association by a vector which is nearly parallel to the extinction vector, could 
the association be only an illusion caused by intervening extinction with Ak ~ 0.25? If so 
we would expect an overall deficit of stars in this direction (due to the extinction of dim 
sources), but there is no such trend. In fact the missing sources are about half as numer- 
ous as the surplus ones. Corroborating thi s point, we note that t he vector separating the 



two populations is slightly steeper than the iNishiyama et al.l (120081 ) extinction law (whereas 



the overpopulation is perfectly consistent with OB stars reddened by that law). Moreover, 
sources which appear underdense on this plot vary on larger angular scales than that of the 
association; their population is mostly a function of galactic latitude, although there are 
somewhat fewer of them toward the association. We therefore conclude that the associa- 
tion does in fact exist, but there may also exist localized extinction which affects the color 
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Fig. 3. — Color-color (top) and color- magnitude (bottom) diagrams of the 2MASS point 
sources. The grey points are all points within the two-degree field, while the red points are 
the sources within the association boundaries. The blue points indicate all sources within 
the color cut. The green diamonds indicate the main sequence at a distance of 9.7 kpc 
with an extinction of Ax = 1.0. The location of an 05V and BOV star are indicated for 
reference. The red line indicates the track of an 05V star at various distances, with the 
numbers indicating the distance in kpc. Similar paths are shown for K3II (black dotted line) 
and BOI (black dashed line) stars. 
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Fig. 4. — Binned color-color (top) and color-magnitude (bottom) diagrams of the 2MASS 
point sources within the candidate association boundaries with the binned background re- 
moved from the diagrams. The cluster color cut is indicated with the black lines. The color 
bars indicate the number of stars in each bin. The white arrow on the color-color diagram 
indicates the extinction vector corresponding to an Ax = 0.5 from the unreddened location 
of O stars. 
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distribution of sources. 

To further investigate the role of patchy extinction, we have obtained additional J 
and H band photometry towards a central 3.9' field in the candidate association using the 
Wide Field IR Camera on the 2.5-metre du Pont Telescope, with limiting magnitudes of 
J = 17.1, and H = 16.3 at a SNR of 3 (J. Radigan, private communication, 2010). We 
extract 906 stars in this field. With this additional photometry, we find stars in the range of 
0.8 < J — H < 0.96, with a continuous density of stars throughout this color range. We also 
find no evidence for any shadowing features which would indicate a strong extinction feature. 
This supports our conclusion that the apparent association is not an extinction feature. 



4. Confusing sources 

In order to identify possible contaminating stars within the color cut, we model the NIR 
colors and magnitudes of template stars of different spectral types and luminosity classes. 
We vary the distances of the model stars to examine which stars meet the cut above the 
2MASS limiting magnitudes. We use reference magnitudes, extinction relationship and the 
average extinction-to-distance ratio as in §3.21 All main sequence OB stars down to B1V 
can fall within the color and magnitude limits. With the 2MASS limiting magnitudes, 
assuming the extinction-to-distance ratio, all subtypes of O stars are visible to 13 kpc. All 
main sequence stellar types later than B1V fall below the magnitude limits. OB giants 
(down to B9III) and K1III to M0III giants can fall within the color and magnitude limits. 
All supergiant stars can fall within the association color cut and magnitude limits. We 
show examples of stars that may fall within the color and magnitude limits in Figu re 



The red-clump giant fea ture, composed of early K giant stars (lLopez-Corredoira et al.ll2002l 



I 



Indebetouw et al.ll2005l ) partially overlaps the association color cut towards the faint end, 



but is not the primary contributor to the background population. 



Using the iRobin et al.l (120041 ) model of the stellar population we produce a synthetic 
catalogue of stars. We simulate the field in the direction and with the size of the Dragonfish 
Association, using the 2MASS limiting magnitudes from £j2j and a diffuse extinction law of 
1.6 Ay /kpc. The model background population underestimates the total number of sources. 
However, we find that the majority of stars within the color cut are early K giant stars 
with masses between 1 and 3 M and distances between 3.0 and 4.5 kpc; we assume these 
represent our contaminating sources. 



5. Conclusions 



We have used the 2MASS point source catalogue to search for the central OB associa- 
tion within the Dragonfish Nebula by analyzing the on-sky stellar density with NIR colors 
consistent with extinguished O stars. We find an overdensity of point sources within the 
bubble at an extinction range of 0.96 < A K < 1.2. Comparing the CCD and CMD of the 
candidate association with the surrounding field, we confirm the presence of the overdensity 
at the chosen color range. We reject an extinction explanation for the overdensity and infer 
that the contaminating star population consists predominantly of K giants 3-4.5 kpc away. 
Our candidate association contains 406 ± 102 members visible in 2MASS with magnitudes 
consistent with O and early-B stars, similar to the ionizing population expected from the 
region's free-free luminosity. We infer from its assumed mass, radius, and maximum age 
that it is likely to be unbound on the whole but may contain a bound core; however further 
investigation is necessary to determine its dynamical state. This is a candidate for the most 
luminous (and consequently the most massive) OB Association in the Galaxy. If confirmed 
by upcoming spectroscopic observations, the Dragonfish Association is an outstanding candi- 
date for multiwavelength examination, as it is a promising laboratory massive star formation, 
stellar dynamics, and feedback processes within the Galactic environment. 

We thank J. Radigan for obtaining NIR data for us. We thank N. Murray, P.G. Mar- 
tin, and R. Breton for the many helpful discussions. This publication makes use of data 
products from the Two Micron All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analysis Center/ California Institute of 
Technology, funded by the National Aeronautics and Space Administration and the National 
Science Foundation. 

REFERENCES 

Benjamin, R. A., et al. 2003, PASP, 115, 953 

Binney, J., & Merrifield, M. 1998, Galactic astronomy, ed. Binney, J. & Merrifield, M. 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 

Conti, P. S., & Crowther, P. A. 2004, MNRAS, 355, 899 

Dutra, C. M., Bica, E., Soares, J., & Barbuy, B. 2003, A&A, 400, 533 

Fall, S. M., Krumholz, M. R., & Matzner, C. D. 2010, ApJ, 710, L142 



Figer, D. F. 2008, in IAU Symposium, Vol. 250, IAU Symposium, ed. F. Bresolin, 
P. A. Crowther, & J. Puis, 247-256 

Gieles, M., & Portegies Zwart, S. F. 2011, MNRAS, 410, L6 

Indebetouw, R., et al. 2005, ApJ, 619, 931 

Kennicutt, Jr., R. C, Edgar, B. K., & Hodge, P. W. 1989, ApJ, 337, 761 

Lopez-Corredoira, M., Cabrera-Lavers, A., Garzon, F., & Hammersley, P. L. 2002, A&A, 
394, 883 

Martins, F., & Plez, B. 2006, A&A, 457, 637 
McKee, C. F., & Williams, J. P. 1997, ApJ, 476, 144 
Mercer, E. P., et al. 2005, ApJ, 635, 560 
Murray, N., & Rahman, M. 2010, ApJ, 709, 424 

Nishiyama, S., Nagata, T., Tamura, M., Kandori, R., Hatano, H., Sato, S., & Sugitani, K. 

2008, ApJ, 680, 1174 

Nishiyama, S., Tamura, M., Hatano, H., Kato, D., Tanabe, T., Sugitani, K., & Nagata, T. 

2009, ApJ, 696, 1407 

Pickles, A. J. 1998, PASP, 110, 863 

Portegies Zwart, S. F., McMillan, S. L. W., & Gieles, M. 2010, ARA&A, 48, 431 
Rahman, M., & Murray, N. 2010, ApJ, 719, 1104 
Roberts, M. S. 1957, PASP, 69, 59 

Robin, A. C, Reyle, C, Derriere, S., & Picaud, S. 2004, A&A, 416, 157 
Salpeter, E. E. 1955, ApJ, 121, 161 
Skrutskie, M. F., et al. 2006, AJ, 131, 1163 



This preprint was prepared with the A AS IATgX macros v5.2. 



